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The scaling properties of selfavoiding polymerized membranes are studied using renormalization 
group methods. The scaling exponent v is calculated for the first time at two loop order, v is found 
to agree with the Gaussian variational estimate for large space dimension d and to be close to the 
Flory estimate for d = 3. 



The statistical properties of polymerized flexible mem- 
branes are interesting and still poorly understood [Q. 
These objects arize either in a collapsed (fractal dimen- 
sion df — 3), a crumpled swollen (2 < df < 3) or a 
flat {df = 2) phase. The physical properties of such 
membranes in three dimensions can be studied by exper- 
iments and computer-simulations. Most of the simula- 
tions find a flat phase due to an induced effective 
curvature term (stiffness) of the membrane . A swollen 
phase with fractal dimension near to the Flory-prediction 
df = 2.5 has been found by exactly balancing curvature 
terms with a long-range repulsive interaction, thus inter- 
preting the swollen phase as a tricritical point j^. The 
experimental results are contradictory. In but not in 
a swollen phase is found. 

An analytical approach inspired from polymer the- 
ory [pl|, which relies on renormalization group and e- 
expansion methods, was initiated in §jl^, where it was 
used to perform calculations at 1-loop order. Its consis- 
tency to all orders in perturbation theory has been es- 
tablished in . In this letter we present the first appli- 
cation of this method to 2-loop calculations and discuss 
the results obtained by this approach. 

The membrane is modeled by a continuum model a 
la Edwards: the embedding of the D-dimensional mem- 
brane in d-dimensional bulk space is described by the 
mappings 6 — > f{x) G R"^. The renormalized Hamil- 
tonian is 



{Vr{x)y+bZ,pi^ 6''{r{x)-r{y)) , (1) 



where b is the dimensionless renormalized coupling con- 
stant, fj, the renormalization momentum scale and e — 
2D — d{2 — D)/2. Physical quantities are calculated 
perturbatively in b. Direct calculations for D = 2 for 
membranes are not possible, since perturbation theory 



is singular; D and d have therefore to be treated as con- 
tinuous variables and an e-expansion must be performed. 
The renormalization factors Z{b, e) and Zi,(b, e) are intro- 
duced in order to subtract the short-distance divergences 
which appear as poles in s at the critical dimension e = 0. 

At order b" one has to to calculate the expectation 
value with respect to the free theory (6 = 0) of n 

bilocal operators, henceforth called dipoles: x' 'y = 

6'^(f(x) — r{y)) integrated over the whole membrane. 
Short-distance divergences occur when dipole end-points 
approach each other. The most important tool to deal 
with these divergences is the multilocal operator prod- 
uct expansion (MOPE) jl^], which describes all possi- 
ble contractions of dipoles to the operators marginal at 
£ = 0. Power counting shows that there are only two 
such operators: the dipole operator and the local opera- 
tor: i(Vr(x))^ = -^x- For instance, the contraction of a 
dipole to a point generates -f with the MOPE coefficient 
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The integral over the relative distance of the two points 
is UV-divergent. As in we use the minimal subtrac- 
tion scheme to subtract these divergences. Introducing 
an IR-cutoff L oc fj,^^ the Feynman-diagram becomes 



all distances 
smaller than L 



(3) 



Our strategy is to keep D fixed and to expand (^ as a 
Laurent series in e, which here starts at e' 
ing by ( I ) p the term of order of ( | 



Denot- 
)l=v the 



counterterms are chosen to have pure poles in e, and 
for instance the first counterterm corresponding to (|2|) is 
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Extending this analysis along the lines of ||T^ leads to the following results (details will be given elsewhere [0). 
second order, the counterterms which render the theory finite are found to be 
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Following , the renormalization group /?-function and 
the anomalous scaling dimension v of r are obtained from 
the variation of the coupling constant and the field with 
respect to the renormalization scale /i, keeping the bare 
couplings fixed. They are written in terms of Z and Z\, 
as 
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l + 5|.lnZ, + f6|;lnZ 



(6) 
(7) 



In order to compute these coefficients, we apply our 
methods of ||l^ , which must be extended in order to deal 
with the double poles in from subdivergences. We 
demonstrate the method using the example oj^the coun- 
terterm associated with the diagram Q = 
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2e 



+ + (8) 



A subdivergence occurs when the single dipole to the 
right of the diagram Q is contracted to a point. Ac- 
cording to the MOPE of [|2|, when this contraction is 
performed /?rs<, the MOPE coefficient factorizes as 




(9) 



This implies that the double pole of (g) is the same as the 
double pole appearing in the product of the counterterms 
associated with the two diagrams on the r.h.s. of d) 
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The factor 1/2 comes from the nested integration [ [l3[ , 
arizing from the fact that the double pole is given by 
the r.h.s. of (||), integrated with the restriction Tl. that 

the distance in ( / is smaller than all the distances in 



-f-; . As a consequence, we can extract the differ- 
ence of the residues of the single poles, ci — ci, by sub- 
tracting to the l.h.s. of (H) a counterterm proportional to 
the r.h.s. of (^), restricted to the domain TZ. In fact, this 
combination has to be calculated, see (0). We get: 
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The careful reader will have remarked that (|1C|) and the last term in (|ll|) are calculated with a slightly different 
regularization prescription. The difference however is of order and thus does not change the result ]T9[ |. ci — ci is 
now extracted by applying to the r.h.s. of (111]) and taking the limit e ^ 0. We then obtain 



ci- ci = ^ lim 

2 e^O 



{largest "1 
diatancc = I- J 



(12) 



This integral is locally finite and can be reduced to an 
integral over five independent distances between pairs of 
points. We evaluate it by applying and extending the 
numerical techniques of |18|] . The main difficulty thereby 
comes from the fact that although convergent, the inte- 
gral ( |l2|) has integrable singularities, which have to be 
removed by suitable variable transformations and map- 
pings between domains of integration. Finally the in- 



tegrand has large variations in some small subdomains 
and can only be integrated by a genuine adaptive Monte 
Carlo integration. For analytical and numerical details 
we refer the reader to |l^ , [l^ . 

The other diagrams appearing in ^ and (||) are 
calculated similarly. Having performed all numerical 
calculations (~ lO'^h CPU on a Workstation), the 
renormalization-group functions f3{b) and iy{b) can be 
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calculated. (3{b) has a non-trivial IR-fixed point for 
b = ba > (i.e. /3(6c) = and P'{bc) > 0). The full 
dimension of the field ly at this critical point, iy(bc) is a 
function of D and e: 



D 



+ z/i(-D)e + i^2(-D)e' + 0(£3) (13) 
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FIG. 1. The functions i^iiD) (dashed line) and V2{D) (solid 
line). The latter is given with the statistical error. 

The coefficients are plotted in figure 1. Naively set- 
ting D — 2 and d = 3, i.e. e = 4 in equation (13), 
yields the wrong result 0. To obtain the correct re- 
sult, one has to redevelop this expansion around any 



point {Do, do 



2-Do 



) on the critical curve (e = 0). 



We use a generalization of the methods introduced in 
pO| . Note that the expansion ( p^ ) is exact in D and 
of order 2 in e, thus can be expanded up to order 2 
in D — Do and e. Given any invertible transformation 
{x,y} = {x{D,e),y{D,e)}, one can express D and e as 
function of x and y and re-expand up to order 2 in a; and 
y. The goal is to find a set of variables, such that the es- 
timate for v depends the least on the choice of the expan- 
sion point on the critical curve and which reproduces well 
the known results in the following cases: D = 1, D = d 
and d ^ oo (see below). The sets {D, Dc{d) = ■^^} and 
{Dc{d),e} have been found to be good choices. The plot 
in figure 2 shows the value of as a function of the expan- 
sion point {D, dc{D) — iD /{2 — D)} on the critical curve. 
The prediction at 1-loop order (dashed line) is essentially 
independent of the expansion point. At 2-loop order the 
estimate starts from the 1-loop result at _D = 2, grows un- 
til it reaches a plateau around D — 1.5 and then grows 
rapidly again. This is a general feature of this kind of 
expansion and can well be studied by applying the same 
method to the Flory-estimate J^Fiory = {2 + D)/{2 + d). 
In this case the plateau becomes flatter if one goes to 
higher orders but does not cover the whole range, i.e. the 
expansion is not convergent for all D. To extract v from 
figure 2, one uses the maximum and the minimum of the 
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FIG. 2. Extrapolation of u via (^3|) in Dc{d) and e for mem- 
branes D — 2 in d — 3. The first order results are given by 
the dashed line, the second order by the solid line. 

plateau. Their mean is an estimate for z/, their difference 
an estimate of the error in this expansion scheme. We 
get I' = 0.74 ±0.02. It must be emphasized that different 
sets of variables yield different values for i^. 

One possibility to further improve the extrapolation is 
to develop vd or i'{d + 2). The first is interesting, as it 
calculates corrections around the estimate predicted by 
a Gaussian variational ansatz |21|, i/var = 2D/d: 



v{b^)d = 2D 



/3(fe)^lnZb(5) 



(14) 



The smaller In(Zb) is, the more accurate the expansion 
becomes. At 1- and 2-loop order h\{Zb) vanishes like 
exp(— const/d) for large d. We argue that this persists to 
any order in perturbation theory. For large d, i^var thus 
becomes exact. 
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FIG. 3. Extrapolation of v via ( |l5| ) in d and e for mem- 
branes D = 2 in d = 3. The first order results are given by 
the dashed line, the second order by the solid line. 
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The second possibility is a systematic expansion 
around the Flory-estimate t'Fiory = {D + 2)/{d + 2): 



v{b^){d + 2) = D + 2 + 



(15) 



b=b^ 



The Flory-estimate is excellent for polymers and one 
hopes that it is also good for membranes. Therefore 
the corrections should be small. Such an extrapolation 
is given in figure 3. The estimate for v is slightly larger 
than that in figure 2. 

The results of a 2-loop extrapolation for v are pre- 
sented in figure 4 for membranes {D = 2) in d-dimensions 
(2 < d < 20). 



is stable towards perturbation by bending rigidity. In- 
deed, it has been argued |3| that for small enough d this 
might destabilize the crumpled phase and explain why 
numerical simulations in d = 3 normally see a fiat phase. 
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FIG. 4. Extrapolation of the 2-loop results in d and e for 
membranes D = 2 in d dimensions, using ( p^ (squares) . The 
solid line is the prediction made by Flory's theory, the dashed 
line by the variational ansatz. 

We see that for d ^ oo the prediction of the Gaus- 
sian variational method becomes exact, as argued above. 
For small d, the prediction made by Flory's argument is 
close to our results. This is a non-trivial result, since the 
membrane case corresponds to e = 4 and in comparison 
with polymers in d = 3, where £ = 1/2, the 2-loop cor- 
rections were expected to be large. In fact we have found 
that they are small when one expands around the critical 
curve e = for an adequate range of £> 1.5 (depending 
slightly on d and on the choice of variables) and a suitable 
choice of extrapolation variables. In this case the 2-loop 
corrections are even smaller that the 1-loop corrections 
and allow for more reliable extrapolations to e = 4. 

In conclusion, we have presented here the first renor- 
malization group calculation at 2-loop order for self- 
avoiding fiexible tethered membranes. In order to im- 
prove these results, one should understand if the plateau 
phenomenon observed at 2-loop order persists to higher 
orders, and one should control the general large order be- 
havior of perturbation theory for this model. Another im- 
portant issue is whether the IR-fixcd point studied here 
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